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We examine the violation of the k± factorization approximation for quark production 
in high energy proton-nucleus collisions. We comment on its implications for the open 
charm and quarkonium production in collider experiments. 

I. Introduction 

Semi-hard processes, where i/i ^ nig^ ^ Aqcd, contribute significantly to particle pro- 
duction in high-energy collider experiments due to the large density of the small-x gluons. 
The k± factorization formalism [ JJ systematically resums corrections of (a^ ln(s/g^))"' 
from gluon branchings in perturbative QCD. In this framework, the particle production 
cross-section is expressed as a convolution of a hard matrix element and unintegrated 
distributions of gluons in the hadrons with definite transverse momentum ki±_ and longi- 
tudinal fraction Xj in each projectile hadron {i=l, 2). 

Multiple-scattering (higher twist) effects become important at small x due to the large 
density of small-x gluons. It is expected to be the origin of the Cronin enhancement and 
p± broadening of hadrons observed in nuclear experiments. It is also relevant for the 
nuclear suppression of quarkonium production. 

The simplest situation for studying the impact of higher twist effects on k± factorization 
is in proton-nucleus (pA) collisions, wherein the proton is dilute and the nucleus is dense. 
The k± factorization formalism was examined in the color glass condensate framework [|2]. 
It is shown that the factorization is recovered when one keeps only the terms that are of 
the lowest order in the charge sources pp^A of the projectiles[^. The cross-sections at the 
leading order in pp, but at all orders in the dense source p^ of the nucleus are obtained 
analytically. Gluon production by the "2-to-l" processes is shown to be /c^-factorizable [ 
mElE] whereas the quark production is generally not [ |Hl E] • 

Here we report the numerical estimates for the k± factorization breaking in quark 
production within the McLerran-Venugopalan (MV) model [ilO|. We briefly discuss open 
charm production and quarkonium suppression in pA collisions in this framework. 
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2. Violation of k± factorization in quark pair production 

The quark pair production cross-section is obtained as[[7]: 

da alN f 6{pj_ + q^- kij_ - k2±) 



d^p^d^q^dypdyg 87i^{N^ - 1) 



+trd 




+m)TgS-m)-f^Ty + h.c.l (Pf'{k,^; k^) 



where the exphcit forms for the Dirac matrices Tgg{ki±, k^) and Tg{ki^) are given in [ 
|7j. Reie (pp{l±) = {n'^Rpg'^ /l^} F.T. (^pp{0)pp{x±)) is the unintegrated gluon distribution 
for the proton, and F.T. denotes the Fourier transformation. One needs, however, three 
nuclear distributions defined as (see Eqs. (42), (43) and (45) in [[?]) 



A 



(/±; fc±) = ^^^7^ F.T. tr (uix^)eU^{y^)t'Uta{0: 



fHh; k±; k'^) ^ ^^^^^ F.T. tr {u{Oru\y^)U{x'^)t'^U^{y'^: 



(2) 



where U and U denote the path-ordered exponentials of the gauge fields in the nucleus 
in the adjoint and fundamental representations, respectively, and describe the multiple 
scatterings of the gluon and the quarks. The average (■ ■ ■ ) is taken over the Gaussian 
distribution of the color charge sources characterized by the saturation scale Q^. 

k± factorization is violated by the transverse structure of the quark pair probed by the 
momentum k±^'^ from the nucleus since each quark from the pair can resolve and interact 
with several gluons from the nucleus. If any of the transverse masses m^^ and of the 
produced quarks is large compared with the typical rescattering scale, Qs, we can neglect 

fcj^ mTqg{ki±,k^!^) and recover the k± factorized formula thanks to the sum rule for 0a 's; 

r j.mm _ r jaw _ ±9,9 

In Fig. 1 we compare the exact result with the /c_l factorized approximation for single 
charm quark production. The breaking is relatively small for the saturation momentum 
(5s=l GeV^, which may be the relevant scale for RHIC at central rapidity. At (5^=15, 
25 GeV^ (corresponding to very forward rapidities in the proton fragmentation region at 
RHIC and LHC) the correction can be as large as 40% at q_L ^ Qs- For the bottom quark 
production the violation is smaller. To assess the model-dependence of our results, we 
compute them now, shown in Fig. 2, with a non-local Gaussian model known to be the 
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exact / kj^-factorized (m = 1.5 GeV) 
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Figure 1. Breaking of factorization in 
single charm quark production. 
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Figure 3. spectrum of the quark pair 
with fixed invariant mass M. 



exact / kj^-factorized (m = 1 .5 GeV) 
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Figure 2. The same as in Fig. 1 but in the 
nonlocal Gauss model. 



exact / kj^-faclorized for pairs (m = 1.5 GeV) 
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Figure 4. Breaking of factorization in 
charm quark pair production. 



asymptotic solution of renormalization equations for x evolution[ llj; non-linear evolution 
effects reduce the magnitude of the violation of k^ factorization. 

In Fig. 3 shown is the total Pj_ distribution of the charm quark pair with the fixed 
invariant masses M=3.1, 4, 8 GeV. In the /c_l factorized approximation (thin curves), 
either quark or antiquark exchanges all the momentum from the nucleus and we see the 
bump structure near Qg, reflecting the gluon distribution of the nucleus. The bump 
is smeared out due to multiple scatterings of both the quark and antiquark in the full 
formula. Integrating over Pj^, we show in Fig. 4, the magnitude of factorization breaking 
in the invariant mass spectrum of the pair. 

3. Phenomenology 

We study the importance of small- a; distributions in D meson production by convo- 
luting the single quark spectrum with an appropriate fragmentation function [ I12j. We 
find, however, the production spectrum is determined not by the quark distribution with 
(iL^Qs-, but largely by the tail part oc of the MV model. Moreover, in order to assess 
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Pair Spectrum, m=1.5GeV 
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Figure 5. dependence of the charm 
pair production. 
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Figure 6. Suppression of low mass pairs 
in pA colhsions. 



the rapidity dependence of open charm production, the x-dependence of the unintegrated 
gluon distributions should be taken into account, which requires going beyond the MV 
model. Our results on open charm production will be reported elsewhere [ [13J. 

The Q^-dependence of the pair spectrum (divided by the charge density is diplayed 
in Fig. 5. At larger M, where the high-density effects are diminished, all curves converge 
to a single one. The multiple scatterings of the pair quarks suppress the yield in the 
low M region. (The overall cross-section is of course enhanced with increasing Q^.) One 
can get an idea about the normal suppression of the quarkonium production in the pA 
collisions, relying on the color evaporation picture. We show the nuclear modification 
ratio, RpA, for the pairs with M less than the open charm threshold 2Md, as a function of 
Q^. The suppression pattern fits the form 1/{Q1)'^ with a ~ 0.42, and not the frequently 
assumed exponential form. One should note here that ~ A^^^ in the MV model. 
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